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We propose a novel conceptual design for a graphene-based quantum engine, driven by a super¬ 
position of mechanical strain and an external magnetic field. Engineering of strain in a nanoscale 
graphene flake creates a gauge field with an associated uniform pseudo-magnetic field. The strain- 
induced pseudo-magnetic field can be combined with a real magnetic field, leading to the emergence 
of discrete relativistic Landau levels within the single-particle picture. The inter-level distance and 
hence their statistical population can be modulated by quasi-statically tuning the magnetic field 
along a sequence of reversible transformations that constitute a quantum mechanical analogue of 
the classical Otto cycle. 

PACS numbers: 05.30.Ch,05.70.-a 


I. INTRODUCTION 

The concept of a quantum heat engine (QHEN) has 
been discussed extensively in the literature [iHlJI as an 
alternative to efficiently recover, on a nanoscale device, 
thermal energy in the form of useful work. In a QHEN, 
in contrast with a classical heat engine, the working sub¬ 
stance exhibits quantum mechanical properties. Interest¬ 
ing examples of this concept are constituted by photosyn¬ 
thesis in plants [l^ , as well as human-designed photocells 
1,0, where the working substance are thermalized pho¬ 
tons. Moreover, it has been recently proposed that if the 
reservoirs are also of quantum mechanical nature, these 
could be prepared into quantum coherent states @,0 or 
into squeezed thermal states [l^ , thus allowing for an en¬ 
hancement of the theoretical photocell efficiency beyond 
the classical Carnot limit 00 0- These examples 
by no means constitute the only possible configurations, 
since a number of different designs based on alternative 
principles have been proposed in the literature, such as 
entangled states in a qubit 0 and quantum mechanical 
versions of the Diesel [I^ and the Otto cycle 0, [H, 0 . 
At the conceptual level, a statistical ensemble of con¬ 
fined single-particle systems can undergo a cycle of re¬ 
versible transformations driven by a generalized external 
field. The driving field can be a mechanical force EHIH 
that modifies the inter-level spacing of the single-particle 
spectrum, thus inducing a sequence of transitions on the 
statistical population of the single-particle states, in close 
analogy to a classical gas confined by a piston. We have 
discussed generalizations of this idea in the context of rel¬ 
ativistic Dirac particles 0] . More recently, we proposed 
01 a magnetically driven QHEN, based on the combined 
effects of a parabolic confining potential, representing a 
semiconductor quantum dot, and an external magnetic 
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field. In the single-particle picture, this configuration 
possesses an exact solution in terms of effective Landau 
levels, which constitute a discrete spectrum, where the 
inter-level distance can be modulated by tuning the ex¬ 
ternal magnetic field 0|. 

In the present work, we propose yet a different alter¬ 
native, through the conceptual design of a QHEN based 
on a strained single-layer graphene flake. It has been 
widely discussed in the literature 0-0| that within 
the single-particle picture provided by the effective Dirac 
theory, mechanical strain in graphene translates into a 
gauge field whose curl plays the role of a pseudo-magnetic 
field. Beyond theoretical calculations, this effect has 
been experimentally observed 0] through the evidence 
of pseudo-relativistic Landau levels in strained graphene, 
in the absence of an actual magnetic field. Moreover, 
via ’’strain-engineering” 0, 0| it is possible to create 
a mechanical strain field that induces a nearly uniform 
pseudo-magnetic field up to Bs ^ 300 T 0. On top of 
this strain field, it is of course possible to impose a real 
uniform magnetic field perpendicular to the plane of the 
sample, with the expected additive effect. In this work, 
we propose a sequence of reversible transformations, in¬ 
volving the quasi-static tuning of the external magnetic 
field, that constitute a quantum mechanical analogue of 
the classical Otto cycle. We calculate the efficiency of 
this cycle, and compare with the Carnot efficiency at the 
same temperatures of the source and sink thermostats. 
In particular, we show that better performance is ob¬ 
tained by combining magnetic and strain fields than by 
adjusting strain solely. 


II. STRAINED GRAPHENE IN AN EXTERNAL 
MAGNETIC FIELD 

The effect of strain-induced pseudo-magnetic fields on 
the electronic properties of graphene has been exten¬ 
sively discussed in the literature 0-0. Moreover, by 
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FIG. 1: (Color online) The deformation vector field 
\x{x\,X 2 ) = {ui,U 2 , 0) defined by Eq.®. 



strain engineering it is possible to generate nearly ho¬ 
mogeneous pseudo-magnetic fields [^, , thus allowing 

for the emergence of relativistic Landau levels, as con¬ 
firmed experimentally The components of the in¬ 

plane pseudo-vector potential As induced by mechanical 
strain are [ 2 l|, [ 2 ^ 

^S.l = {uil — U22) , ^S,2 = (—2^12) . ( 1 ) 

la la 

Here, the local displacement vector is defined as u = 
[ui,U 2 ,z), with Ui the in-plane components, and z the 
out-of-plane component, while /3 = 91nt/cHna is the rel¬ 
ative change of the hopping parameter 6t/t with respect 
to the dilation of the lattice constant 6a/a [ 2 ^. The 
strain tensor [uij] is defined by 

^ + dizdjz). ( 2 ) 

It was recently shown [23, that the continuum 
expansion of the tight-binding Hamiltonian for strained 
graphene in the vicinity of a single Dirac point (valley), 
up to 0{ufj), is given by the expression 


be expressed in the so-called ’’laboratory” frame, that is 
with all the fields described with respect to the instanta¬ 
neous atomic positions in the deformed crystal configura¬ 
tion yi = Xi-\-Ui{x). This statement leads to the inclusion 
of the last contribution into Eq. (|4]) for the velocity tensor 
[ 2 ^, given by Uij = ^ {diUj + djUi). On top of the gauge 
field A 5 , a ’’geometric” vector field T^ appears, whose ori¬ 
gin can be understood within the covariant picture [13] 
as the pseudo-spin connection for fermions propagating 
in a curved space. When expressed in the ’’laboratory” 
frame [l^, it becomes 

~ ”"4 (5) 

The pseudo-magnetic field generated by the 
strain-induced gauge potential is given by 
Bs = 63 ( 9 iHs ,2 — d 2 As,i). The last term in the 
Hamiltonian Eq.Q represents a pseudo-Zeeman cou¬ 
pling between the pseudo-spin degree of freedom and 
the strain pseudo-magnetic field Bs [ 2 ^ . 

Let us consider the deformation field (see FiglT]) 

ui = -2usxiX2, U 2 = -us{x\- x\), z = 0, ( 6 ) 


where us is a constant characterizing the magnitude of 
the in-plane displacement, in association through Eq. (HD 
with the vector potential 


As = 2^^ {-X2, xi, 0). 


This in turn implies a pseudo-magnetic field 
Bs = 4- 63 , 


(7) 

( 8 ) 


which is constant in magnitude and points along the out- 
of-plane direction 63 . Let us now discuss what is the 
effect of this particular deformation on the local Fermi 
velocity. First, it is straightforward to check that Ukk = 
V • u = 0, and that since z = 0, we have Uij = Uij. The 
local Fermi velocity tensor is then directly calculated by 
substituting these relations into Eq.(|4]) 


v^j = VF 


2 ^b 1 — XF6k 


(9) 


+VFaiTi - jBsas] ip{x). (3) 

Here, vf 10® m s~^ is the Fermi velocity for the undis¬ 
torted graphene lattice. The presence of strain induces a 
local position dependence on the velocity [U, [ 23 ] , which 
then becomes a tensor 

Xijix^ — XF {/2Uij SijUkk') T ■ (4) 

De Juan et al. have shown that in order to describe 
correctly the effects of strain, the Hamiltonian should 


where in the last step we used that for graphene P = 
dlnt/dlna ^ 2 [U, [^. Therefore, for the deformation 
field defined by Eq. (| 6 D , the velocity is a constant diagonal 
tensor Vij = VFSij. 

Finally, let us analyze the corresponding expression for 
the field T^ dehned by Eq.fO. After Eq.dHD we have Ujj = 
0, and hence diUjj = 0. Moreover, it is straightforward to 
verify from Eq.(l 6 D that djUij = djUij = 0. Therefore, for 
the deformation field defined by Eq.([nD we have T^ = 0 , 
and the Dirac Hamiltonian Eq.(I3]) reduces to 

H= d^x{x) {vF^ ■ [—i'^ + As] — pBsas} tjj{x), 

( 10 ) 
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with a- = ((Ti,(T 2 ) and As given by Eq.©. 

The Hamiltonian in Ea. ifTU)) describes the physics in 
the vicinity of a single Dirac point (valley). The inclu¬ 
sion of the pseudo-magnetic field terms breaks the valley 
degeneracy since, as opposed to a real magnetic field, the 
strain-induced gauge potential has opposite signs in the 
vicinity of each Dirac point [25j |. Therefore, we define 
the Hamiltonian describing both valleys by the higher¬ 
dimensional structure 

H = J d^xi’Hx) 

Here, = ± represents each of the valleys, corre¬ 
sponding to the two inequivalent points the 

first Brillouin zone, respectively. By adopting the stan¬ 
dard convention for the components of the spinor field at 
each sublattice (A, B) 

( 12 ) 

the effective Hamiltonian at each valley, according to 
Eg. lfTTI) is given by 



The Hamiltonian Eg. dTSl) determines a system of two 
decoupled Dirac equations, one for each electronic spin 
component s = {±}, 

{^vf^ ■ [—zV -I- -I- A] — syHl — ^^Bsas) ipi^\x) 

= (17) 

The pseudo-spinor eigenstates of the Hamiltonian in 
Eq-dIZl) are with (see Ap¬ 

pendix A for details) 



Here, the coefficients q;„ and f3n are as defined in Ap¬ 
pendix A (see Eg. llAlSI) ). The corresponding energy 
eigenvalues are 


Eis 


-I- A| — sjB, 

- sjB, 


n > 0, 
n = 0. 


(19) 


= J (fxip^^'^'^{x){^VF^-[—i^ + ^A.s] 

-^^Bsd3}ip^^\x). (13) 

Let us now consider the combined effect of a uniform 
magnetic field B = e^B in addition to the strain-induced 
pseudo-magnetic field defined defined by Eq.®. Choos¬ 
ing for convenience the gauge 

A = y (-a;2,a;i,0), (14) 

we have that the effective Hamiltonian at each valley Aj, 
including the electronic spin degree of freedom becomes 

= J d^a;'I'^^^’^(a;) (1 0 O') ■ [—zV-I-^As-I-A] 

-^iBs (1 ® da) - yH (da 0 1)} 'I'^^^(a;). (15) 

Here, we have defined the spinor 

#(€) [x) = ^ 0 ^i^^(16) 

where is the two-component pseudo-spinor at valley 
A'f arising from the bipartite graphene lattice, whereas 
the electronic spinors Xf = (l^O)^ and Xi = ( 0 : 1 )^ are 
eigenvectors of da with eigenvalues s = {±}. The last 
term in the Hamiltonian Eg. (1151) is the Zeeman interac¬ 
tion, characterized by a coupling constant 7 = g/ZB/ 2 , 
with g ~ 1.8 for graphene [^. Notice that the Zeeman 
coupling only involves the real magnetic field B, since the 
strain-induced pseudo-magnetic field B5 does not inter¬ 
act with the electronic spin. 


Here, A = ± is the band index, while = VF\/‘2e\B^\/h 
is the effective frequency, expressed in terms of the ef¬ 
fective ’’total” magnetic field B^ = B + ^Bs that re¬ 
sults from the combination of the strain-induced pseudo- 
magnetic field Bs and the real magnetic held B at each 
valley K^. We have also dehned the dimensionless pa¬ 
rameters Aj = at each valley. The two 

quantum numbers (n, m) correspond to the quantiza¬ 
tion of the orbital Landau level n > 0, and the guiding- 
center (i.e. the center of the classical cyclotronic orbit) 
coordinate m > 0, respectively [ 2 ^. The energy lev¬ 
els described by Eg. (|TOl) are degenerate in the guiding- 
center quantum number m, with the same degeneracy 
factor ^ for each Landau level and each val¬ 
ley. On the other hand, A| = is the number 

of magnetic flux quanta (jjn = /z/ 2 e piercing the area A 
of the graphene flake [^, with $ 5 ^ = B^A the ’’to¬ 
tal” hux. As previously discussed, only the magnetic 
field B couples to the electronic spin s = {±}, as seen 
in the Zeeman term in Eg. (1191) . which is proportional 
to 7 = ( 7 ^b/ 2 . Regarding the pseudo Zeeman term, 
when Bs is expressed in Tesla, we estimate [2l|, [H, [ 2 ^ 
7 = = 9.788 X 10-^eVT-i, for a = 1.42 A the 

carbon-carbon bond length, and V = 6 eVA [11,133. 
Remarkably, 7 ~ 1.7^b is on the order of magnitude of 
the Bohr magneton. 

III. MAGNETIC QUANTUM ENGINE 

As the ’’working substance” for our quantum en¬ 
gine, let us consider a statistical ensemble of repli¬ 
cas of a single-particle system [S Hi consisting on 
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FIG. 2: (Color online) The cycle is pictorially represented in 
the entropy (S) versus external magnetic field (B) coordinates. 
The cycle is composed by two iso-entropic trajectories, and 
two trajectories at constant external magnetic field. The cold 
reservoir is at Ti = Tc, whereas the hot reservoir is at Ts = 
Th. 

a single electron in the conduction band (A = +) 
of to the graphene flake described by Ea. dTSl) . This 
can in principle be achieved by charging an otherwise 
neutral graphene flake with a positive gate potential. 
Each replica [l^, may be in any of the different 
eigenstates of the Hamiltonian Eq. (II3- The single¬ 
particle system is then in a statistically mixed quan¬ 
tum state [3l[, described by the density matrix op¬ 
erator p = Yjn,rn,s,iPri,rn,s,i{B)\^l}^n}rnAB)){'4’n}rnAB)\, 
with \ipn)m,siB)) a spinor eigenstate of the single-particle 
Hamiltonian Eg. dTTl) for a given magnetic field intensity 
B and pseudo-magnetic field Bs- The indexes (n, m, s, 
enumerate the eigenstates of the Hamiltonian defined by 
Eg. lfT^ . with Landau level n, at valley ^ and electronic 
spin component s, for A = -I-. This density matrix opera¬ 
tor is stationary, since in the absence of an external per¬ 
turbation [3l| ihdtp = [H, p] = 0. Here, the coefficient 
0 < Pn,m,s,^{B) < 1 represents the probability for the 
system, within the statistical ensemble, to be in the par¬ 
ticular state \il^n}m,siB)). Therefore, the {pn,m,s,^{B)} 
satisfy the normalization condition 

Trp = ^ Pn,m,s,({B) = 1. (20) 

In the context of Quantum Statistical Mechanics, en- 
tropy is defined according to von Neumann as 

S = — fcsTrplnp. Since in the energy eigenbasis the equi¬ 
librium density matrix operator is diagonal, the entropy 
reduces to the explicit expression 

S{B) = -ks ^ Pn,m,s,^{B)\n{pn^ra,s,i{B)) . (21) 

In our notation, we emphasize the explicit dependence 
of the energy eigenstates {|'!/’i^L,s(H))}, as well as the 


probability coefficients {pn,m,s,^iB)}, on the intensity of 
the external magnetic field B. 

The ensemble-average energy E = {H) of the quantum 
single-particle system is 

E = Tr(piL) = ^ p^,m,e,dB)Eis{B) 

= Y.P^^sAB)EiAB), (22) 

n,s,^ 

where we introduced the coefficients 

Pn,s,^ = ^ ^ Pn,Tn,s,^ 
m 

in order to take notational advantage of the degeneracy 
in the spectrum with respect to the quantum number m. 

The statistical ensemble just described can be submit¬ 
ted to an arbitrary quasi-static process, either by modu¬ 
lating the magnetic field intensity or by exchanging en¬ 
ergy with a reservoir. Along such process, the ensemble- 
average energy will change accordingly [ 2 , , 

dE = Pn,sAB)dEls{B) + dPn,sAB)EijB) 

n,s,^ n,s,^ 

= ' (^ 4 ) 

The first term in Ea. dMl) represents the energy change 
due to an iso-entropic process, whereas the second term 
represents the change due to a process where the energy 
spectrum remains rigid. These two terms are in cor¬ 
respondence with the macroscopic notions of work and 
heat, respectively. Therefore, Eo. d^^ represents a mi¬ 
croscopic version of the first law of thermodynamics for 
the statistical ensemble of single-particle systems [mis, 
dE = SWj^s + SQ^^s- Here, according to Eo. (l24l) the 
work along the process connecting states with magnetic 
fields By — 7 ^ Bs is 



On the other hand, the heat exchanged by the system 
with the environment while modifying its temperature 
from Ty —>• Ts will be 

For the statistical ensemble just defined, let us con¬ 
sider a cycle by devising a sequence of quasi-static tra¬ 
jectories in Hilbert’s space, as depicted in Figl^] Initially, 
the single-particle system, while submitted to an exter¬ 
nal magnetic field of intensity Bi, is brought into thermal 
equilibrium with a macroscopic thermostat at tempera¬ 
ture Ti = Tq. In equilibrium, von Neumann’s entropy 
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achieves a maximum, given by the Boltzmann distribu¬ 
tion (/3 = {kBT)~^) 

Pn,s,dB) = [Z{B,p)]-^ (27) 

where the normalization factor is defined by the partition 
function (see Appendix B for details) 

Z(i?,/3)= (28) 


The ensemble-average energy for the statistical distribu¬ 
tion in Eg. d^Tl) is given by the expression (see Ea. (IB6l) in 
Appendix B) 

= <»> 

It also follows from the definition Eq. d^Tl) . along with 
Ea. (l?fl) and the normalization condition Eg. (1201) that the 
entropy can be expressed in terms of the partition func¬ 
tion (see Eo. dBTp in Appendix B) 


S{l3,B)/kB = (3E{l3,B)+\nZ{B,l3) 

= -/3(^)^+lnZ(i?,/3). (30) 

The system performs work along the iso-entropic tra¬ 
jectory 1 —?► 2, according to Eo. d^ . 

/ B 2 / f)F'\ 

=E{T2,B2)-E{Tc,BiX31) 

and along the iso-entropic trajectory 3 —4, 

W3^4 = dB =E{T4,B4 )-E{Th,B2). 

(32) 


A physical interpretation of the work performed by the 
engine is obtained by considering the statistical mechan¬ 
ical definition of the ensemble-average magnetization, 
that is M = — {dE/dB)g. Hence, the work terms de¬ 
fined in Eo. d^Tl) and Eg. (15^ can also be interpreted as 
W = -f MdB. 

Along the constant magnetic field trajectories 2 —?► 3 
and 4 —?► 1, the system exchanges heat with the reservoirs. 
The heat absorbed by the system from the ’’hot” reservoir 
at Ts = Th is 



= E{Th,B 2 )-E{T 2 ,B 2 ). (33) 


Similarly, the heat released by the system to the ” cold” 
reservoir at Tq is 


Qc 



E{Tc,B4)-E{T4,B4). 


(34) 



FIG. 3: (Color online) The efficiency of the cycle, as a function 
of the compression ratio r{B 2 ), for the case 7 = 0 (red, dash- 
dot line) compared with the case 7 = 1.7/rs (blue, solid line). 
Here Bi = 4T, Bs = 20 T, and the temperatures at the 
reservoirs Th = 100 K, Tc = 30 K, respectively. 


The efficiency of the engine is then given by the ex¬ 
pression 


V = 


Wi ^2 + W3^4 


Qh 


= 1 - 

Qc 


Qh 


(35) 


The intermediate temperatures T 2 and T 4 must be deter¬ 
mined numerically from the condition that connects the 
initial and final states along each iso-entropic trajectory 
(see Figl2]) 


S{BuTc) = S{B2,T2), 

S{B 2 ,Th) = SiB4,T4). (36) 


Eor given values of the initial magnetic field Bi, the 
strain pseudo-magnetic field Bs, and the reservoir tem¬ 
peratures Tc and Th, the efficiency is a function of the 
magnetic field i? 2 . We choose to parametrize this depen¬ 
dency by defining the ratio 


r{B 2 )=lBjlB„ (37) 

where Ib = min{Z+,l_} is a characteristic confinement 
length for the semi-classical cyclotronic orbit, defined as 
the minimum Landau radius among the two inequivalent 
valleys. In Fig. [3] we plot the result of our numerical cal¬ 
culation of the efficiency, as a function of the magnetic 
field expressed in terms of the ratio r{B 2 ). In this partic¬ 
ular example, we have chosen Bi = 4 T, Bs = 20 T and 
the temperatures Tq = 30 K and Th = 100 K at the cold 
and hot reservoirs, respectively. In Fig. [3l we also com¬ 
pare the effect of the pseudo-Zeeman term, by calculating 
the efficiency when setting 7 = 1.7pB and 7 = 0, respec¬ 
tively. It is evident from comparison of both curves that 
the pseudo-Zeeman effect produces a relative enhance¬ 
ment of the efficiency as compared to the case when this 
term is absent. 
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We find that the numerical solution for the efficiency 
as a function of the compression ratio can be accurately 
represented by the parametric form 

r?=l-[r(B2 )]-“. (38) 

Here, the exponent a depends on the temperatures Tq 
and Tff, as well as on the strain field Bs- In particular, 
for the choice of parameters represented in FiglHl we find 
that a = 6.88 for 7 = 1 . 7 ^ b , whereas a = 4.2 for 7 = 0. 
Remarkably, this parametric expression is analogous to 
the well known formula for the efficiency of the Otto cycle 
that works with a classical ideal gas, with r instead of 
the volumetric ratio that applies to the classical case. 
An even closer analogy between both cases can be put 
forward by recognizing that r(i?2) = Ibi/Ib 2 > 1 can 
be literally interpreted as a ” compression ratio” between 
the effective Landau radii, which in practice defines a 
characteristic confinement length for the semi-classical 
cyclotronic orbit associated to each Landau level. 

Numerical solutions for the efficiency only exist up 
to a maximum compression ratio, which in the exam¬ 
ple displayed in Fig. 0 is r"*®® = 1.19 for 7 = 1.7/r_B, 
whereas = 1.33 for 7 = 0. At this point, the effi¬ 
ciency attains its maximum value, that exactly matches 
the Carnot efficiency for the same temperatures in the 
thermostats, i.e. = r]c = I — Tc/Th = 0.7. 

More generally, using the parametric form Ea. (l38L one 
concludes that 

r<(^y (39) 


IV. CONCLUDING REMARKS 

We believe that for a practical implementation of the 
cycle proposed in this work, it should be more conve¬ 
nient to apply a combination of a real magnetic field 
and strain, than with strain solely. Experimentally, it 
is easier to impose a dynamic modulation of the applied 
magnetic field intensity, rather than to modulate strain 
over a nanoscale sample. In our proposed scheme, only 
a static strain pattern is imposed. On the other hand, 
using only a real magnetic field in the absence of mechan¬ 
ical strain, presents the disadvantage of requiring very 
high magnetic field intensities. This problem is avoided 
by the strain induced pseudo-magnetic field, which ex¬ 
perimentally can achieve up to Bs ~ 300 T [2^, and 
thus provides a high baseline for the total effective field 
\B± I = \B±Bs \ that allows for the emergence of relativis¬ 
tic Landau levels, which as we propose can be modulated 
with a real magnetic field B of moderate intensity. 
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Appendix A 

In this appendix, we provide details of the solution for 
the single-particle spectrum of the Hamiltonian Eq. (II3, 
following the operator method presented in [2^. Let us 
first define the canonical momentum operators at each 
valley AT^, 

n« = -ihv 4- e (A -f ^As). (Al) 

The single particle Hamiltonian at each valley for 
^ = ± is therefore expressed by 

— ^^Bso's — s'yBl. (A2) 

The eigenvalue problem can then be formulated in the 
form 


{^VF^ ■ n« - 'y^Bsas) (A3) 

where we have defined e| = A| -|- s^B. Let us now de¬ 
fine the Landau radius associated to the total effective 
magnetic field B^ = B + ^Bs at each valley, 


I 


« = 



(A4) 


It is straightforward to verify that the components of the 
canonical momentum operator satisfy the commutation 
relation 

[ni, n«] =-i^ sgn{B^), (A5) 

which allows us to introduce the creation and annihila- 


tion operators 



a = 

isgn(H^)ng), 


= 

(ni + *sgn(H^)ng), 

(A6) 

satisfying the canonical commutation relation 



[a, o’^] = I 

(A7) 


as a consequence of Ea. (|A6l) . In explicit matrix notation, 
the eigenvalue problem then reduces to (for sgn{B^) > 0) 




r-A^ 

a 

( \Un) ^ 

- 

\nu) \ 

at 


V y 

^n,s \ 

\Vu) ) 


(AS) 


with Aj = . The matrix structure for 

sgn(Hj) < 0 is obtained under the similarity transfor¬ 
mation and hence it shares exactly the 

same spectrum. After Ea. (IA8l) . we have the eigensystem 
of equations 

^hn^a\vn) = |w„), (A9) 

^nn^d^\un) = \vn). (AlO) 
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For the ladder operators defined by Eas. (IA6l) and (ITtI) . 
the combination a^a = fi is the number operator. Com¬ 
bining Ea. dAQI) and Ea. dAlOl) . we have 


a^a\vn) 




(All) 


Clearly then, |u„) = |n) is an eigenstate of the number 
operator, a^a\vn) = a^a|n) = n|u„). Therefore, from 
Ea. dAll|) we obtain 


+ A|, n > 0, (A-12) 

where the sign A = ± arising from the two roots of the 
quadratic equation determines the valence and conduc¬ 
tion bands, respectively. The first component |u„) of the 
spinor, for n > 0, is determined from Ea. dAQI) . using 
a\vn) = y/n\n — 1), as follows 


\Un) = 




Cn.s 


y/n 


Aj -I- Y^n -I- A| 


a\vn) 

'n - 1) (A13) 


The corresponding (normalized) spinor eigenvectors 
are thus given by 




where we have defined the coefficients 

^/n 


I3n = 


-I- A| -I- A^Aj -I- A|^ 
A^ -I- A^ Y^n -I- A| 

^71 -|- A| -|- A^Aj y^ti -f A|^ 


1/2 ' 


1/2- 


(A15) 


The case n = 0 must be analyzed separately. From 
Ea. dAQI) . we conclude that |u„) oc a|u„) oc |n — 1). 
In particular, since |u„=o) = |0)i then it follows that 
= 0. Applying this condition to Ea. dAlOl) for 
n = 0, one immediately concludes 


eL = (A16) 

The corresponding spinor eigenvector for the Landau 
level n = 0 is thus given by 

IV'Eo.s) = ( |S) ) ’ ^ = 0- 


So far, we have only considered the quantization of 
the orbital degrees of freedom. However, an intrinsic 
(very large) degeneracy in the eigenvalues described by 
Ea. dA18|) emerges when we take into account the quan¬ 
tization of the guiding-center of the orbit, i.e. the center 
of the classical cyclotron orbit. For that purpose, one 
decomposes the position operator r = (x,y) = R -I- p 
into the guiding-center coordinate R = (A, Y) and the 
cyclotron orbital variable p — {px, Py)- The later is per¬ 
pendicular to the electron’s velocity, and hence it is ex¬ 
pressed in terms of the canonical momentum components 


Px 




(A19) 


In consequence, one has the commutation relations that 
follow after Ea. dASI) 


[Px, Py] 


[n|,n|] 

ieBif 


-ill sgn(R^). 


(A20) 


The commutation between the components of the posi¬ 
tion operator [x, y] = 0 implies the commutation relation 

[X,Y] = - [px,Py] = IZ|sgn(H^). (A21) 

Therefore, we are entitled to define a second pair of cre¬ 
ation and annihilation operators 

b = -^{X + isgn{B^)Y), 

= ^^{X-zsgn{B^)Y) (A22) 


that satisfy the canonical commutation relation [6, (jf] = 
1 , thus defining a pair of ladder operators, characterized 
by the eigenvalue equation 

b^b\m) = m|m). (A23) 


Therefore, the full solution of the eigenvalue problem 
stated in Ea. (IA3l) is obtained by the tensor product of 
the spinor eigenstates defined in Eas. (IAI4l) . (IA17I) and 
the eigenstates |m) associated to the guiding-center co¬ 
ordinates 


\^P 


(«) 

X,n,m,s 


) 


IV’S,.) ® \^) 


1 f an\n - I,m) 

^2 V Pn\n,m) 


n > 0 


l '0 


(0 X 

n=0,m,s/ 


l^i= 0 .s) ® k) 


0 

|0,m) 


n = 0. 

(A24) 


The spectrum is then obtained by the shift ^ ^ = 

A,n,s - siB, 

^ I Xhfl^^n -f A| - syB, n > 0, 

I — sjB, n = 0. 


Here, the coefficients an and /3„ are as defined in 
Eg. dAISp . The uncertainty relation between the guiding- 
center coordinates that follows from the commutation re¬ 
lation Ea. (IA21l) is [2^ 


AAAF = 27r/|. (A25) 
































Just as in the case of phase-space coordinates, this uncer¬ 
tainty relation defines the minimal area of an elementary 
surface cell associated to a given eigenstate. Therefore, 
the degeneracy ^ of each Landau level, corresponding 
to the maximum value of m allowed, is given by the num¬ 
ber of elementary cells contained in the total surface A 
of the sample 


The degeneracy of each relativistic Landau level n > 0 
with respect to the center of mass quantum number 
m > 0 corresponds to the number of effective mag¬ 
netic flux quanta piercing the area A of the sample, 
^n,s = -^1 = Under these considerations, the 

partition function factors into the form 


= Ni = ——— 
AXAY 


A 

27rl| 


4>o 


(A26) 


Here, = B^A is the flux of the ’’total” effective 
magnetic field piercing the area of the sample, whereas 
(j)o = h/e is the universal magnetic flux quantum. Notice 
that in general, the degeneracy factors for a given exter¬ 
nal magnetic field B and strain pseudo-magnetic field Bs 
will be different at each valley K^. 


Appendix B 


^ = 2cosh(/37i?) (iV+Z+ + N-Z.) . 


Here, we have defined 


_ g-j/3KnjA^ 


OO _ 

n—1 



(B3) 


(B4) 


The single-particle spectrum corresponds to the rel¬ 
ativistic Landau levels for the effective magnetic field 
B^ = B + ^Bs, involving the strain Bs aX each valley 
and the real magnetic field B contributions 


^ I fin^^n + A^- s-fB, n>0, 

I — sjB^ n = 0, 

for = ^Bshr^fl^^. Here, we consider only particle¬ 
like solutions (A > 0), assuming that the graphene flake 
has been charged with a single electron in the conduction 
band. The partition function that describes the statisti¬ 
cal ensemble of single-particle systems is given by 


—1 

OO ^^n,s 

^ = EEE E ' 

s—dz n—0 m—0 

OO 

= EEE<^ 

s=± n—0 






(B2) 


For the low temperatures of interest, and at relatively 
large magnetic fields, one finds that at both valleys 
^ 1. For example, consider T ^ 100 K and 
B^ ~ 10 T, which are of the order of magnitude of the 
numerical examples discussed in the main text. For this 
case, one readily finds (3hn^ ^ 13.3. Given this fact, the 
exponential terms which add up to the infinite series are 
very small, and hence the sum can be well approximated 
within less than 2% of relative error by retaining only the 
first term, 




(B5) 


Using this expression, explicit analytical formulas are ob¬ 
tained for the ensemble-average energy E 


Ei/3,B) = - 


InZ 

/ B 


= —^B tanh(/37i?) 




(B6) 


as well as for the entropy S/ks = In Z — 

= /3U(/3, B) + In [2 cosh(/37H)] (B7) 


E^l 

5 =± 


_j_ g-/3fif2{.yT+A| 


-I- In 
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